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Abstract 

We calculate the flux of 7-rays emitted from primordial black holes (PBHs) 
which are formed by a "blue" power-law spectrum of density fluctuations in 
the early universe. Gamma-ray emission from such PBHs may contribute 
significantly to the observed extragalactic diffuse 7-ray background (DGB). 



Using the observed DGB flux from the imaging Compton Telescope (COMP- 
TEL) and the Energetic Gamma Ray Experiment Telescope (EGRET) as the 
upper limit of 7-ray flux from PBHs, we derive the upper limit on the spectral 
index n of the density fluctuations. The range of initial PBH masses which 
can contribute to the DGB is 2 x 10 13 g — 5 x 10 14 g, corresponding to a cosmic 
reheating temperature of 7 x 10 7 GeV — 4 x 10 s GeV. In this range, we find the 
upper limit to be n < 1.23 — 1.25. This limit is stronger than those derived 
from the energy density in PBHs or PBH relics and matches the value of n 
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required to explain the cosmic microwave background anisotropy. 
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I. INTRODUCTION 



Distinct from black holes which form by recent processes like stellar collapses, black holes 
formed by mechanisms in the early universe can exist. Such black holes, named primordial 
black holes (PBHs), produce many interesting consequences in the early universe and can 
also be sources of present astrophysical events. The simplest mechanism for PBH formation 
is the density fluctuations in the early universe [|I|]. Overdense regions which strongly deviate 
from the background universe can evolve into black holes when the overdense regions enter 
the cosmological horizon. The resulting PBH mass is about the horizon mass at entry. 
Thus, PBH masses can be as small as about the Planck mass Mpi ~ 2 x 10~ 5 g or as large 
as 1O 7 M . In the latter case, they can bound the mass of a typical galaxy after decoupling 
@. PBHs surviving today, or their massive relics, can be sources of dark matter. 

Particle emission from black holes due to Hawking evaporation enlarges the role of 
PBHs. Interactions of the emitted particles with the matter in the universe can affect 
numerous early universe phenomena, such as nucleosynthesis |4j], baryogenesis ||, cosmic 
microwave background radiation (CMBR) distortion ||, entropy production 0, diffuse 7-ray 
background (DGB) [||[8|-|nj, and so on. PBHs with initial mass ~ 5 x 10 g are presently at 



the final stage of their evaporation |TTJ and may emit enormous amounts of energy. In some 
GUT-scale theories, the relics of PBHs created with initial mass < 5x 10 14 g and expired by 
today can constitute the dark matter [12|. 

From the effects of PBHs mentioned above, upper limits on Qbh, the fraction of the 



critical energy density of the universe which can be in PBHs, are found |13| . Constraints on 
the spectral index of the density fluctuations have been derived from these energy density 
limits [|14|,|T5| and it has been shown that PBHs in general give the strongest upper bounds 



on the spectral index |16 



The upper limit of Qbr in PBHs with masses about 5 x 10 14 g was also calculated from 
the PBH 7-ray emission [§,§. There exists a homogeneous and isotropic 7-ray background 



in the universe whose origin is known to be extragalactic [17,Tq]. Earlier authors postulated 
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a contribution to the DGB from the PBH 7-ray emission |5|f|[|. MacGibbon and Carr 
have recently updated the limit on the PBH density using the DGB measurement of the 
Energetic Gamma Ray Experiment Telescope (EGRET) [B| and found that Qbh ~ (5.1 ± 
1.3) x io- 9 /i- 195±0 - 15 g. Here h is the Hubble parameter in units of lOOkms^Mpc" 1 . 
Similar values are deduced from the PBH antiproton and 7-ray emission if PBHs cluster 
along with cold dark matter in galactic halos |2ll , |22l . The limits on Qbr in this mass range 



assume, though, that the density fluctuations have a scale-invariant Harrison-Zel'dovich 
spectrum with spectral index n = 1. Therefore they cannot be converted into an upper 
limit on the spectral index. Also, in these approaches the fluctuation amplitude was not 
explicitly normalized. Instead, a parameter related to the PBH density was introduced and 
varied to match the PBH 7-ray flux to the observed DGB. However, if one normalizes the 
fluctuation amplitude on the scales for PBH formation to that detected by COBE on much 
larger scales from the CMBR anisotropy amplitude S ~ 1.9 x 10~ 5 [23], it is impossible to 



form a significant number of PBHs with a continuous Harrison-Zel'dovich spectrum. With 
the normalization to the CMBR anisotropy, significant PBH abundance is possible only if 
the density fluctuations have an n > 1 ( "blue" ) spectrum. A blue spectrum with a constant 
spectral index is a valid assumption, for example, in the hybrid inflationary scenario |TJJ . 



For these reasons, we will re-examine the 7-ray flux from PBHs formed by an n > 1 
spectrum and find the upper limit on the spectral index, using the recent DGB measurements 
of the imaging Compton Telescope (COMPTEL) @] and the EGRET @ on board the 



Compton Gamma Ray Observatory (CGRO) and normalizing the fluctuation amplitude to 



that on CMBR anisotropy scales. While previous authors [|14i|16l have estimated limits on 
n from the 7-ray emission of PBHs formed by a blue fluctuation spectrum, they have not 
performed the explicit calculation of the PBH 7-ray emission and matched it to the DGB, 



analogous to the approach of MacGibbon and Carr for n = 1 PJZUf. They also did not 
include the effect of quark and gluon production which dominates the emission above black 
hole temperatures of about lOOMeV . 

In Sec. II, we review the PBH formation and PBH mass spectrum, correcting errata 



in Ref. |15| . The black hole evaporation and QCD fragmentation effects on the particle 
emission are discussed in Sec. Ill and the calculation of the 7-ray flux from PBHs is given 
in Sec. IV. The recent DGB observations are reviewed in Sec. V. Our detailed calculation 
of the upper limit on n is presented in Sec. VI. The paper closes with some concluding 
remarks in Sec. VII. 

II. PBHS AND THEIR MASS SPECTRUM 

We address PBH formation in a universe with a hard equation of state, that is p = 
with < 7 < 1. Studies of the evolution of a spherical overdense region whose initial radius 
R is greater than the particle horizon show that for the region to collapse to a black hole, 



the initial density contrast of the region, Si, should satisfy the following condition [26 



Here a and (3 are constants of the order of «/t~, Mj is the mass contained in the region of 
radius R at the time t% when the fluctuation develops, and Mm is the horizon mass at t{. 
The lower bound comes from the requirement that the radius of the region at its maximum 
expansion should be larger than the Jean's length at that epoch. The upper bound comes 
from the requirement that the overdense region should not be disconnected from the universe. 
Since Mj oc R 3 and R oc k~ l , we will use Mj, R and k interchangeably to represent the initial 
mass, size or comoving wavelength. The PBHs form when the overdense region enters into 
the horizon. The resulting PBH mass, Mem, is approximately the horizon mass at that time 
£h an d is given by 

M B ffi^7 3/2 M Hi ^ . (2) 



Mbhi is related to Mi via [2B| 



M Bffi ^ 7 37/{1+37) M i {1+7)/(1+37) M^ /{1+37) . (3) 
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As the universe expands, larger PBHs are formed, so that PBHs with masses less than Mem 
coexist in the universe at time £h- 

We assume the density fluctuations to be Gaussian. Recently it was claimed that a 



non-Gaussian nature to the fluctuations may affect the PBH formation [P?]. This is model 



dependent and does not much alter the upper limit on n ]T5| . If one surveys the universe 



with a window having size R, the smoothed density field 5r(x) is defined by 

<5fl(x) = Jd 3 y5(x + y)W R (y) , (4) 

where d>(x) = (p(x) — pb)/pb, Pb is the background energy density of the universe, and Wr(x) 
is the smoothing window function of size R. The dispersion ctr, the standard deviation of 
the density contrast of the regions with R, is given by 



where Vw ~ -R 3 denotes the effective volume filtered by Wr, and 5k and are the Fourier 
transforms of 5(x) and W^x), respectively. For Gaussian fluctuations the probability that 
the region of size R has density contrast in the range (5 + d5, 5) is 



1 r d 3 k 



1 / 5 
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P(M i ,5)d8=-=^exp[-- T )d5 . (6) 



Thus, the probability that the region with Mj collapses to a black hole is 

Pbh(M 4 ) = [ A p(M t ,5)d5 , (7) 



with 



The above quantity Pbu(M) has been interpreted as the ratio of the density in PBHs to 
the density of the universe. This is not strictly so because regions larger than M, also 
contribute. As in our previous work [[DJ, we proceed by omitting these contributions and 
find the number density of black holes produced by the collapse of regions with mass between 
Mi and M { + dM { to be @ 
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TlBuiM^dMi = ~ 



Pi 12 B da R 



Mi\ na 2 R dMi 



exp 



B 



2a 



(9) 



Here pi = 3/(32nGt 2 ) is the background energy density of the universe at t{. Unlike Ref. 



15 1, we will not convert the mass spectrum into a function of unsmoothed quantities but 



retain the use of the smoothed quantities. This is to avoid factors being dropped in the 
conversion. 

Since fluctuations grow as a 2 (t) in the radiation-dominated era, where a(t) denotes the 
cosmic expansion factor, the dispersion corresponding to <tr becomes at horizon crossing 



0"H 



OR 



(10) 



whence nBu(Mi) can be written as 



nBE(Mi)dMi 



TtM 7 



1 <9cr H 2 CTjj 1 



a\ dMi 3 Mi 



exp 



(id 



Applying Eq. (3) with 7 = 1/3, the mass spectrum can be expressed as a function of the 
initial PBH mass Mbh« and we have 

-3 r 1 <9cr H 



^BH(M B ffi)rfM B ffi 



■^Ip^f^M^ 



0", 



exp 



7 

2^ 



BHi 



(12) 



.o&dM Bm M Bm _ 

In general, tensor perturbations (gravitational waves) are also produced in inflationary 
scenarios and contribute to the CMBR anisotropy. However, inclusion of the tensor per- 
turbations does not significantly affect the PBH mass spectrum if the fraction in tensor 



perturbations is not dominant ||14|| . Therefore we assume that the anisotropy is only due to 
the scalar fluctuation. 

The four-year results of the COBE experiment resolve the horizon crossing amplitude at 



present to be |23j 



S = 1.91 x l(T 5 exp(1.01(l -n)) 



(13) 



with n — 1.2 ± 0.3. This implies a smoothed amplitude today a of 9.5 x 10 5 with a slight 
dependence on n |TB[ . We denote the mass contained at ij in the region whose comoving scale 



corresponds to the present horizon scale by Mq. Under the power-law spectrum assumption, 
cr H oc k^ 1 ^ 2 and so 



(Ml 

VMo 



l-n 
6 



(14) 



(Note that the spectral index n of Ref. [26| is equivalent to (n + 3)/6 in this work.) Mq is not 
the present horizon mass as it was incorrectly taken to be in Ref. That misidentification 
introduced further errors when converting o"h into a function of M^m and made the results 



of Ref. [15] far weaker than those of Ref. [16]. Our revised results are given in Sec. VI. 
From Eq. (3), oh can be represented as a function of M Bm 

( M Bm y 



cr H = (To 



I: 



(15) 



v-^BHiO ' 

withp= (1 — n)/4 in the radiation-dominated era (7 = 1/3) and p= (1 — n)/6 in the matter- 
dominated era (7 = 0). Minor discrepancies, which arise in Eq. (15) at the transition into 
the matter-dominated era, lead to less than a 1% change in the constraint on n [IB]. 

The PBH initial mass spectrum under the n > 1 power-law spectrum assumption is thus 
described by 



n BH {M BHi )dM BUi 



n + 3 



7 

2al 



BHi 



(16) 



III. PARTICLE EMISSION FROM BLACK HOLES 

Due to the Hawking effect || , a rotating charged black hole emits particles at a rate, 

6XP l n/2n J + ( - 1) J ' (17) 



dN s r s 



duodt In 

per degree of particle freedom. Here, k, Q, and $ are the surface gravity, angular velocity 
and electric potential, respectively, s is the particle spin, I is the axial quantum number or 
angular momentum and q is the particle charge. The absorption probability for the emitted 
species, T s , is in general a function of u, Q, k, together with the internal degrees of freedom 
and rest mass of the emitted particle. Taking the above emission rate, it has been shown 
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that — > before most of the black hole evaporates |28j and that a black hole with mass 
< 1O 6 M discharges faster than it evaporates |]29"fl . Hence, it is natural to regard PBHs as 
Schwarzschild black holes. 

At high energies, T s oc u 2 for massless or relativistic particles and the emission spectrum 
mimics the radiation from a black body of temperature Tbh = k/2t[. Noting that the surface 
gravity is k = 1/4GMbh when Q = <3> = 0, the temperature of a black hole can be defined 
as 

T BH = * ^ 1.06(^f) 'GeV . (18) 
87tGM B h V10 13 g7 v ; 

At low energies, T s does not simply scale as u 2 but depends on other quantities mentioned 
above. The form of T s has been explored both analytically and numerically (see references 
in [p0[] ). Extensive numerical studies of the direct particle emission from black holes were 
done by Page ||31|| . Hawking emission can be thought of as a process by which a black 
hole emits particles with approximately a blackbody radiation spectrum once the black hole 
temperature exceeds the rest mass of the particle. Thus, black holes with masses larger than 
10 17 g, corresponding to Tbh — O.lMeV, emit only massless particles. As the black hole mass 
decreases, massive particles will be emitted. 

In the conventional viewpoint, it is natural to assume that elementary particles like 
quarks and gluons, rather than composite hadrons, are directly emitted from black holes 
once the emission energy exceeds the QCD confinement scale, Aqcd- In this picture, pions 
are only directly emitted from black holes in the energy range between lOOMeV and Aqcd 
and are produced by quark and gluon decay above Aq CD . Taking into account the number 



of the emitted species, the mass loss rate of a black hole can be written as [1 I 



dt 



-5.34 x lO 25 0(M BH )M B ^g sec" 1 , (19) 



where 0(Mbh), a function of the number of directly emitted species, is normalized to unity 
for Mbh ^ 10 17 g. Black holes with masses 5 x 10 14 g -C Mbh "C 10 17 g emit e ± , neutrinos and 
photons and have initially 0(Mbh) = 1-569. If black holes can emit three lepton families, 
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six quark flavors, the photon and direct pions, then 0(M B h) ~ 13.9. Including the emission 
of weak gauge and higgs bosons, 0(M BH ) ~ 15.4 for T BH < lOOGeV and M BH ^ 10 11 g. At 
higher energies or in non-standard models like supersymmetry or superstrings, 0(M BB ) rnay 
be greater but in general remains less than 100. 

Integrating Eq.(19), the black hole lifetime, r evap , is found to be [|T1[ 

r cvap ~ 1.2 x 10 3 — 

0(M BH ) 

= 6.24 x lO" 27 M^ H 0(M BH )~ 1 sec . (20) 

The jet-like fragmentation and hadronization of the quarks and gluons evaporated above 
Aq CD drastically change the observable spectrum of emitted particles. The evaporated 
quarks and gluons fragment into further quarks and gluons which then compose themselves 
into hadrons on distances greater than Aq^ d in the jet frame. These particles further decay 
into the astrophysically stable particles - photons, neutrinos, electrons and protons and their 
antiparticles. MacGibbon and Webber have shown that this picture is analogous to the decay 
of quark and gluon jets in e + e~ accelerator events and calculate the instantaneous flux of 
particles for 0.2GeV < T BH £ lOOGeV by convolving Eq. (17) with the HERWIG QCD jet 
code P5| . Their results differ strongly from those of previous works which omitted QCD 



emissions and particle decays. They find that the black hole emission at these temperatures 
is dominated by the jet fragmentation products. In the case of photons, the primary peak 
in the black hole emission is due to the decay of jet-produced tt° and occurs around 67MeV. 
The position of the photon peak does not shift significantly with the black hole temperature. 
In contrast, the photons directly emitted by the black hole, not resulting from jet decay, 
appear at ~ 5T BB with fluxes 4 to 5 orders less than the flux at the jet-dominated peak. 
The jet-produced photons were omitted in previous estimates of the 7-ray emission from 



PBHs formed by n > 1 density fluctuations [fn,IB]|. In this paper, we fit the instantaneous 



7-ray fluxes from T BB = 0.2 — lOOGeV black holes which are shown in Fig. 4 of Ref. [25 
and derived including quark and gluon emission. We then use these results to calculate the 
7-ray flux from PBHs. 
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IV. GAMMA-RAYS FROM PBHS 



To calculate the flux of 7-rays from PBHs, it is important to know how many PBHs have 
existed in the universe. This can be determined from the PBH initial mass spectrum, which 
in turn depends on the fluctuation amplitude and spectral index, and £j, the time when 
the fluctuations develop. Previous works considered only the case in which the fluctuations 
are described by the Harrison-Zel'dovich spectrum and hid the effect of fluctuation ampli- 
tude and ti [f^SHIOH . However, it is impossible to form a significant number of PBHs with 
the Harrison-Zel'dovich spectrum if the fluctuation amplitude on PBH formation scales 
is normalized to the amplitude found by the COBE experiment on much larger scales. 
With normalization to the CMBR anisotropy, substantial PBH formation is possible only 
if the fluctuations increase on small scales. This occurs if the fluctuations satisfy an n > 1 
power-law spectrum. We will assume that the fluctuations follow a power-law spectrum and 
calculate the 7-ray flux from PBHs with the initial mass spectrum given by Eq. (16). 

Since the number density of PBHs decreases as -R -3 , the number density of PBHs at the 
time £1 > £h is 



\ -3 rM Bm 

^bh(^i)- (-ft) / n BU (M Bm )dM BHi , (21) 
\rLjJ JM*(t\) 



where M*(ti), the initial mass of a PBH whose lifetime is t%, is given by 

<f>(M*(ti)) ( h mV3 



(22) 



.6.24 x 10- 27 Vlsec 

We denote the instantaneous 7-ray flux from a black hole with mass Mbh as / 7 (Mbh,^)- 
At £1, PBHs with initial mass M Bm have evaporated down to a mass 

M evap ~ (M 3 Hi - 1.6 x lO 26 0(M Bffi )ti) 1/3 (23) 

and are emitting photons with flux / 7 (M evap , to). The angular frequency at emission, lu, is 
redshifted by the expansion of the universe to a present angular frequency ujq of 

w = |%, • (24) 
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Thus the total flux per unit solid angle at today, to, of 7-rays emitted from PBHs is 

JJ- = ^ L U) U) * /m.<„, -)»bh(Mbh,)<;M bh , (25) 

where £ m j n is the earliest time after inflation at which PBHs form. 

Photons emitted by the black holes may interact with ambient matter in the universe 
via many processes and lose energy or be cut off during propagation. If a photon effectively 
interacts r times during flight, the photon flux is attenuated by a factor of e~ T . Therefore, 
the actual photon flux reaching Earth at the present time is 

dJ 1 r f o /.Rn\/-Ri\~ 3 r M Bm 
duj Q Att Jt min \RiJ \RiJ JM,{t x ) 

The number of interactions r, known as the optical depth, depends on the energy and the 
time or redshift z at emission. A detailed treatment on the optical depth is given in Sec. 
VI. 



V. EXTRAGALACTIC DIFFUSE GAMMA-RAY BACKGROUND 



Since its first discovery in the 0.1 — 2MeV range by detectors on the lunar probes [p2]| , the 



homogeneous and isotropic diffuse 7-ray flux, whose origin is extragalactic, has been observed 
in numerous satellite and balloon-borne experiments. The SAS-2 satellite provided the first 
clear evidence for the existence of an extragalactic 7-ray background between 30 — 150MeV 
pr§|| . Several Apollo and balloon-borne experiments also saw evidence of a bump in the few 
MeV range in excess of the extrapolated X-ray continuum |33 . 



A number of models were proposed to explain the early measurements of the extragalactic 
spectrum. Active galaxies, which can be observable sources of discrete extragalactic 7-ray 
emission when they are located close to our Galaxy, are believed to contribute at least in 



part P^}| . Another model which has been considered is matter-antimatter annihilation at 
the boundaries of superclusters PSfl . In this model, the MeV bump is attributed to the 
redshifted peak of 7r° decays at 67MeV. PBH 7-ray emission has also been proposed as a 
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contributor to the DGB flux. In some exotic models, PBH emission may additionally explain 
the MeV bump |J. However, none of the scenarios for the DGB production, by themselves, 
is sufficient to explain the measured flux and spectrum. 

With higher sensitivity and wide-field of view, the detectors on the CGRO enlarge the 
detection range and gather important data on the DGB. The DGB flux measured by the 
COMPTEL [[24J at 1 — 30MeV is now compatible with power-law extrapolations of the 
measured flux at lower and higher energies. The results below about 9MeV are preliminary 
but the 2 — 9MeV flux is far less than previously measured and no MeV bump is seen in this 
region at the levels reported previously. This weakens the need to explain an MeV feature. 



We parameterize the preliminary COMPTEL results [24| in the range 0.8 — 30MeV by 



the best-fit power law function given in Ref. [36 



dJ 



due 



obs 



/ u \ -2.38 

6.40 x 10~ 3 ( y^^J [cm 2 s sr MeV]" 1 (COMPTEL) . (27) 



In the range 30MeV - lOOGeV, the EGRET experiment finds the DGB flux to be well 
described by the single power-law function ]T9"| |, 



dJ 



duJr 



= k 

obs 



U 



451MeV/ 



) [cm 2 s sr MeV]" 1 (EGRET) (28) 



with k = (7.32 ± 0.34) x 10" 9 and a = 2.10 ± 0.03. No large scale spatial anisotropy 
or deviations in the energy spectrum is discernible in the extragalactic component above 
30MeV. The observed flux above lOMeV, and possibly up to lOOGeV, may be explained by 



unresolved blazers [fL9|| . Below lOMeV, the measurements still have large uncertainties and 
the exact nature of the emission is not well understood. 



VI. CONSTRAINTS ON THE SPECTRAL INDEX 

We now derive the constraints on the spectral index from the condition that the PBH 
7-ray flux should not be larger than the observed DGB flux. With the normalization of 
the fluctuation amplitude to the CMBR anisotropy, PBH formation is limited to the epoch 
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when the fluctuation arises. This time is related to the reheating temperature in inflationary 
models by |J7| 

t , m = , mg: ^^y^ . (29) 

Here g* ~ 100 counts the degrees of freedom of the constituents in the early universe. The 
minimum initial PBH mass corresponding to Trh is 



1 ,/,._ /T, 



-2 



z RH ^^7 3/2 M P1 ^j . (30) 

PBHs created before the onset of reheating will be diluted to an insignificant density during 
inflation. Since we are considering the case n > 1, the resulting PBH initial mass spectrum 
has a very narrow mass range. Thus we will make the approximation that the photons are 
solely emitted by PBHs whose initial mass is Mrh. 

The photon interactions in the matter-dominated era, relevant to the DGB observations, 
are Compton scattering, pair production, photo-ionization, and photon-photon interactions. 
Via these processes and cosmo logical redshift, the energy of the emitted photons is degraded. 
Zdziarski and Svensson have studied the attenuation of 7-ray flux at cosmological distance 
p8| . They found that the maximum redshift from which photons can be detected today 
peaks at z max ~ 700 and for present energies lMeV < u < lGeV . All photons emitted 
at higher redshifts are cut off by interactions (i.e, t(uq,z) > 1) and do not reach Earth. 
This means that PBHs which completely evaporated before z max ~ 700 can not contribute 
to the present DGB. From Eq. (20), this corresponds to a minimum detectable initial PBH 
mass of about 2 x 10 13 g and a reheating temperature of Trh — 4 x 10 8 GeV. Noting that 
M*(to) — 5 x 10 14 g, the range of PBHs which can contribute to the observed DGB is then 
2 x 10 13 g — 5 x 10 14 g and the corresponding reheating temperature range is 7 x 10 7 GeV — 
4 x 10 8 GeV. Outside the range lMeV ^ c^o ~ lGeV, the maximum redshift for a given 
energy, z max (u ), is less than 700 and depends on the details of the interactions. We take 



z max (uJo) from Ref. ||38|| . Only photons for which r < 1 are here included and we regard 
these photons as being free from attenuation. The integrated 7-ray flux from PBHs does 
not depend significantly, though, on the details of the optical depth. 
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We now proceed to calculate the integrated PBH 7-ray flux, Eq. (26), using the in- 
stantaneous emission from individual black holes, / 7 , obtained by fitting the simulations 



of Ref. |[25|| . Our results are shown in Fig. 1. It can readily be seen that the 7-ray flux 
from PBHs can not fully explain the observed DGB flux although the PBH emission may 
contribute significantly to the observed DGB flux around 10 — lOOMeV. In the n > 1 case, 
the PBH flux arises from the lifetime emission of PBHs with initial mass Mrh, whereas the 
Harrison-Zel'dovich spectrum produces a broad range of initial PBH masses. In both cases, 
the PBH flux falls off as roue hly luq 3 above lOOMeV g,|,l. This high energy tail mainly 
comes from the lifetime direct photon emission in the most recent evaporation epoch |TT] . 



At low energies, where the flux is strongly determined by QCD jet fragmentation, the flux 
spectrum for n > 1 does not scale as Uq 1 , as for n = 1, but instead flattens out due to the 
narrow PBH initial mass range. In addition, the turnover in the spectrum occurs at lower 
energies as the reheating temperature increases and Mrh decreases. This is because the 
turnover corresponds to the redshifted peak emission of an Mrh black hole emitting at its 
initial temperature. 

From the constraint that the PBH 7-ray flux can not be larger than the observed DGB 
flux, we derive the upper limits of the spectral index in the range 7 x 10 7 GeV < T RH ~ 
4 x 10 8 GeV (Fig. 2). The upper limit on n is 



n 



< 1.23 - 1.25 . (31) 



Even though photons emitted in models with higher Trh suffer more interactions and larger 
redshifts, the exponential dependence of the PBH mass spectrum, Eq. (16), implies that the 
number density of PBHs strongly increases with Trh for a given spectral index. Thus, the 
upper limit on n decreases as the reheating temperature grows. That is, the constraint on 
n becomes stronger as the reheating temperature increases. These values are similar to the 
upper limits obtained from the deuterium destruction constraint in the lower mass range 
10 9 g < M BH < 10 13 g 1T§. 



Also shown in Fig. 2 are the weaker limits on n derived from the maximum allowable 
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energy density in PBHs. We plot the upper limits found from the requirement that the 
PBH energy density does not overdose the universe at any epoch, Qbr < 1 (Case I, the 
dashed line), and the similar restriction on any present relic density in PBHs which did not 
evaporate completely but left residual masses of about the Planck mass, fi re ii c < 1 (Case 
II, the dotted line). The latter constraint strengthens somewhat if the relic mass is greater 
than Mpi [|l4l , |T5l . Constraint Case I applies regardless of whether PBHs evolve into massive 
relics. The new upper limits on n from the energy densities are much tighter than those 



of Ref . |1J,|15|] and decrease as Trh grows. In Ref . [16j , upper limits on n were found from 
the condition that the present PBH density and any relic density satisfy ^bho < 1 an d 
^reiic < 1) respectively. There it was shown that the constraint ^bho < 1 is weaker than 
the relic constraint, Case II. However, when we now extend that condition to require that 
the PBH density fraction at any epoch does not overdose the universe, we find that Case 
I is stronger than Case II if Trh ~ 10 13 GeV. The energy density in PBHs or PBH relics 
give weaker constraints than the DGB because the upper limit on ^bh from the PBH 7- 
ray flux is far less than 1. Because of the vast difference between the scales on which the 
CMBR anisotropy occurs and the scales on which PBHs form, the limit on the spectral index 
obtained from PBH emission is highly insensitive to the true value of the CMBR anisotropy 
or Q BU - 



VII. CONCLUSIONS 

In this paper, we calculate the 7-ray flux from PBHs formed by density fluctuations in the 
early universe and compare it with the observed extragalactic DGB flux. Previous works 
considered the case in which the density fluctuations have an n = 1 Harrison-Zel'dovich 
spectrum and did not explicitly normalize the fluctuation amplitude. If the fluctuation 
amplitude on the scales of PBH formation is normalized to that on large scales deduced 
from the COBE observations of the CMBR anisotropy, PBHs can not form in cosmologically 
significant numbers from a Harrison-Zel'dovich spectrum. Thus, we describe the fluctuations 



16 



by an n > 1 power-law spectrum and find the upper limit on n from the condition that 
the 7-ray flux from PBHs should not be larger than the DGB flux. The smallness of the 
fluctuation amplitude from the CMBR anisotropy limits PBH formation to the time when 
the fluctuations develop. In inflationary models, this time is related to the reheating time. 
To model the 7-ray emission, we fit previous simulations of the emission from individual 
black holes which included QCD fragmentation and decays. Due to the interactions of 7- 
rays with the background matter in the universe, only PBHs surviving later than z max < 700 
can contribute to the DGB flux observed today. The initial mass range of PBHs relevant 
to the observed DGB is then 2 x 10 13 g ^ Mbh ~ 5 x 10 14 with a corresponding cosmic 
reheating temperature between 7 x 10 7 GeV < T RH ^4x 10 8 GeV. We find the resulting 
upper limit on n in this range to be n < 1.23 — 1.25. Our constraint on n is stronger than 
those obtained by requiring that the energy density in PBHs does not overdose the universe 
at any epoch (^bh < 1) and that found by requiring that any present PBH relic density 
similarly does not overdose the universe (f2 re i ic < 1). The upper limit on Qbh implied by the 
PBH 7-ray flux is far less than 1. If the fluctuation amplitude is constrained by the CMBR 
anisotropy and PBH emission, the upper limit on n is fine-tuned and highly insensitive to 
the precise upper limit on ^bh or the precision in the CMBR measurement. 



Recently, Niemeyer and Jedamzik [[39(] have argued, supported by preliminary numerical 
simulations, that sub-horizon mass PBHs may form in considerable numbers at any forma- 
tion epoch. For Gaussian fluctuations, they deduce that the PBH initial mass distribution 
at a given formation time peaks at about 0.6 times the horizon mass and extends from much 
smaller masses up to the horizon mass. Such a distribution would have an effect on our 
limits similar to raising Trh- 

We also note that it has been proposed that the emitted quarks and gluons from a black 
hole may interact and form a photosphere around the black hole above black hole tempera- 
tures of a few GeV ||40j| . This scenario, however, remains controversial. In the photosphere 
model, the flux would be more concentrated around lOOMeV than if the quarks and gluons 



directly fragment into hadrons ||41|| . The high energy tail from the PBH distribution would 
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also scale asw 4 , not uj 3 . Photosphere formation may somewhat weaken the constraint on 
n but such changes would be small due to the fine-tuned nature of the constraint. 
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FIGURES 

FIG. 1. The integrated 7-ray flux from PBHs, ^ in units of (cm 2 s sr MeV) -1 , for (a) 
Trh = 7 x 10 7 GeV corresponding to Mbh = 5 x 10 14 g, and (b) Trh = 10 8 GeV correspond- 
ing to M BH = 2 x 10 14 g. The bold lines are the observed flux, ^|obs, from the COMPTEL 
(lMeV <cj < 30MeV) and the EGRET (30MeV < lo < lOOGeV). 

FIG. 2. The upper limits on the spectral index. The solid line between 7 x 10 7 GeV — 4 x 10 8 GeV 
is obtained from the condition that the PBH 7-ray flux should not exceed the observed DGB flux. 
The dashed line is obtained from the condition that Obh < 1 throughout the history of the universe 
(Case I). The dotted line is obtained from the condition that f2 re ii c < 1 (Case II). 
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